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ABSTRACT 

It is an indisputable fact that solar magnetic fields are force-free in the corona, 
where force free fields means that current and magnetic fields are parallel and there is 
no Lorentz force in the fields. While the force-free extent of photospheric magnetic 
fields remains open. In this paper, the statistical results about it is given. The vec- 
tor magnetograms (namely, B x , B y and B z in heliocentric coordinates) are employed, 
which are deduced and calibrated from Stokes spectra, observed by Solar Magnetic 
Field Telescope (SMFT) at Huairou Solar Observing Station (HSOS) are used. We 
study and calibrated 925 magnetograms calibrated by two sets of calibration coeffi- 
cients, that indicate the relations between magnetic fields and the strength of Stokes 
spectrum and can be calculated either theoretically or empirically. The statistical re- 
sults show that the majority of active region magnetic fields are not consistent with the 
force-free model. 

Subject headings: Magnetic field, Photosphere, Corona 
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1. Introduction 



Magnetic fields dominate most solar activites such as filaments eruption, flares and coronal 
mass ejections (C MEs). AH these phenomena are energetic events due to explosive release of 



magnetic energy (Krall et al. 



2001 



Priest & Forbes 



2003 



1982; 



Shibata et al. 



Nindos 



1995 



Tsuneta 



1996 



Wan 



1994 



Bond et al. 



20121) . Thus the solar magnetic field is the key for us 



to understand the nature of solar activities. Since the measurable magnetic sensitive lin es are 
around the photosph ere, the reliable measurements of magnetic fields are nearly there (|Stenflo 



1973 



Harvey 



1977) . However the understanding of the magnetic field in the chromosphere and 



corona remains difficu 



(Garv&Hurford 



as a force-free (Alv 



1994 



t due to both int rinsic physical difficulties and observational limitations 



; Lin et al. 



2004|) . Generally, magnetic fields in the corona are regarded 



1989), because the plasma J3 (ratio of plasma pressure to magnetic pressure 



over there) is much less tha n unity. But it is contro versial in the photosphere, because two kinds 



of pressure are comparable dDemoulin et al. 



1997|). For the low-/? corona where the plasma is 



tenuous (J3 «; 1), magnetic field satisfies the following force-free equations: 



V x B = a(r)B, 



V • B = 0. 



(1) 



(2) 



They imply that there is no Lorentz force in action and a is constant along magnetic field lines 
(B-Va = 0). 

At present, the magnetic field extrapolation with a force-free assumption is a major method 
to study the solar magnetic fields of active regions. The coronal fields can be reconstructed 



from a physical model (namely, force fr ee mode 



field is taken as a boundary condition (IWu et al. 



Sakurai 



1981 



He & Wang 



Yan & Sakurai 



2008: 



Liu et al. 



2000; 



) in which the observed photospheric 



1990; 



Wheatland et al 



Mikic et al. 



2000; 



1994 



Wiegelmann 



Amari et al 



2004; 



magnetic 



1997 



Song et al 



2006 



201 II) . It means that coronal magnetic fields are considered to be 
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force-free, while at the boundary it is connected to the photospheric magnetic fields observed. The 
force-free extent of the p hotospheric magnetic field then becomes an important subject to study. 



Wiegelmann et al. 



(2006) proposed a preprocessing procedure to make a minor regulation within 
the allowable errors, so that the observed magnetic fields are biased to a force-free, which further 
indicates the s tudy o f force-free extent is significant and necessary for the field extrapolation. 



Metcalf et al. 



(I1995|) calculated dependence of the net Lorentz force in the photosphere and low 
chromosphere on the height using Mees Solar Observatory magnetograms and concluded that 
the magnetic fields are not force free in the p hotosphere, but becomes force-free roughly 400 



km above the photosphere. 



Moon et al 



(|2002|) studying the force-free extent in the photosphere 



using 12 vector magnetograms of three active regions, realize d that the photospheric magnetic 



fields are not very far from the force-free case. 



Liu et al. 



(|20 1 IP) tentatively applied the force-free 



extrapolation to reconstruct the magnetic fields above the quiet region and checked the force-free 
extent of this quiet region based on the high spatial resolution vecto r magnetogra ms observed by 



Tiwari 



(|2012|) found sunspot 



the Solar Optical Telescope/Spectro-Polarimeter on board Hinode. 
magnetic fields are not so far from the force-free. We, in this paper, conduct a the statistical 
research to make use of the vector magnetograms observed by SMFT at HSOS from 1988 to 2001 
in order to verify the force-freeness of the photospheric magnetic field. 

The paper is organized as follows: firstly, the description of observations and data reduction 
is arranged in Section|2l then, the results are shown in Section[3l finally, in section |4] presents the 
short discussions and conclusions. 



2. Observations and Data Reduction 

The observational data used here are 925 vector magnetograms corresponding to 925 
active regions observed from 1988 to 2001 by olar Magnetic Field Telescope (SMFT) installed 
at Huairou Solar Observing Station (HSOS), located at north shore of Huairou reservoir. 
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Magnetograms associated with the corresponding to the active region, which is the nearest to disk 
center is chosen and calculated. So only one magnetogram is available for one active region. 
SMFT consists of 35 cm refractor with a vacuum tu be, birefringent filter, CCD camera including 



1986|) . The birefringent filter is 



an image processing system operated by computer (|Ai & Hul 

o o 

tunable, working either at the photosphere line Fe I A 5324.19 A, with a 0.150 A bandpass, or at 
the chromosphere line, HB, with a 0. 1 25 A bandpass. The line of Fe I A 5324. 19 A (Lande factor g 
= 1.5) formed around the solar photosphere, is used for photospheric magnetic field observations. 
The bandpass of the birefringent filter is about 0.15 A for Fe I A 5324.19 A line. The center 
wavelength of the filter can normally be shifted -0.075 A relative to center of Fe I A 5324.19 
A to measure o f the longitudin al magnetic field and then the line center is applied to measure the 



transverse one (Ai & Hu 



1986). Vector magnetograms are reconstructed from four narrow-band 
images of Stokes parameters (/, Q, U and V). V is the difference of the left and right circularly 
polarized images, Q and U are the differences between two orthogonal linearly polarized images 
for different azimuthal directions, / is the intensity derived from either the sum of two circularly 
polarized images is the line-of-sight field measurements or of two linearly polarized images in the 
transverse field measurement. When /, Q, U and V are measured, the corresponding white light 
images are simultaneously obtained, which are employed to compensate for the time differences 
during the measurements of /, Q, U and V. The sequence of obtaining Stokes images is as follows: 
First acquired the V/I image ; next the Q/I images, then the U/I image. The time required to obtain 
a set of Stokes images was about 45 second. Each image is associated with 256 integrated frames. 
To reconstruct the vector magnetograms, the linear relation is necessary between the magnetic 



field and the Stokes t 


parameters /, Q, U and V, 


Jjefferies et al. 


1989; 


Jefferies & Mickev 


1991) 



B L = C L V, 
B T = C T (Q 2 + U 2 ) l/ \ 



(3) 
(4) 
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8 = arctan(-), 
B± 

1 U 
(p = -arctan(-), 



(5) 
(6) 



where B L and B T are the line-of-sight and transverse component of the photospheric field, 
respectively. 6 is the inclination between the vector magnetic field and the direction normal to 
the solar surface and cf> is the field azimuth. C L and C T are the calibration coefficients for the 
longitudinal and transverse magnetic fields, respectively. Both theoretical and empi rical methods 



are used to calibrate vector magnetograms (|Wang. Ai & Deng 



1996 



Ai et al. 



1982). So that two 



sets of calibration coefficients are avai 



(ISu & Zhang 



2004; 



Wang. Ai & Deng 



able: the first set C L and C T are 8381 G and 6790 G 



1996Q . respectively, obtained by theoretical calibration; 
in the second set Cl a nd Ct are 10000 G and 9730 G, respectively, which are deduced through 
empirical calibrations (|Wang et al.l ll996). Fa raday rotat i on an d magneto optical effects may 



affect the value of measured magnetic fields. 



Bao et al. 



(200 0j) analyzed th e Faraday rotation 



Gao et al. 



(120081) pro posed a 



way o f 



Zhangl fcOOOh 



in the HSOS magnetograms which contributes about 12 c 

the statistical removal of Faraday rotation in vector magnetograms from HSOS. 
found that the statistical mean azimuth error of the transverse field is 12.8° caused by magneto 
optical effects. In addition, there may be other causes of the data uncertainty. For example, 
saturation effects, filling factors and stray light. After routine data processing of HSOS data, the 
spatial resolution of observational data is actually 2 arcsec/pixel x 2 arcse/pixel and 3<x noise 
levels of vector magnetograms are 20 G and 150 G for longitudinal and transverse components, 
respectively. 



Wang 



he acute angle method is employed to resolve 180° am biguity (|Wang. Zhang & Ai 



1997 



Wang. Yan & Sakurai 



2001; 



Metcalf et al. 



19941 : 



2006), in which the observed field is 



compared to the extrapolated potential field in the photosphere. The orientation of the observed 
transverse component is chosen by requiring -90° < A0 < 90°, where A0 = 9 -8 e is the angle 
difference between the observed and extrapolated transverse components. At last, the components 
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of the vector magnetic field (B x = B T cos(f>, B y = B T sin(f) (B T = JB^. + Bp and B z = B L cos6) are 
calculated in local heliocentric coordinates. To minimize the projection effects, the requirement 
that the horizontal width of an active region is less than 20 degree is added for each magnetogram. 



3. Results 



In classical electromagnetic theory, Lorentz force can be written as the divergence of the 
Maxwell stress: 



F = 



(B V)B V(B B) 



An 8tt 

Under the assumption that the magnetic field above the plane z = (namely, the photosphere) 
falls off fast enough as going to infinity, the net Lorent z force in the infinite half-space z > is 



(7) 



just the Maxwe 



Alv 



1984; 



Low 



1 stre ss integrated over the plane z = (IChandrasekhar 



19611: 



Molodenskv 



1974 



1985). Then, the components of the net Lorentz force at the plane z = can be 



expressed by the surface integrals as follows: 

F x = --^ J B x B z dxdy, 
1 r 

F y = B y B,dxdy, 

An J 



= ~f (B 2 z -B 2 x -B 2 y )dxdy. 



(8) 

(9) 
(10) 



The necessary conditions of a force-free field are that all three components are much less 



than F p (ILow 



19851) . where F p is a characteristic magnitude of the total Lorentz force and can be 



written as: 



(Eft + B 2 X + B 2 y )dxdy, 



(11) 



The values of F x /F p , F y /F p , and F-JF p then provide a measure of the force-free extent at the 
boundary plane z = (the photosphere). 
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Following Metcalf et all (119951) and Moon etal.1 (120021) . F x /F p , F y /F p , and FJF P are 
utilized to check the force-free extent of the selected photospheric magnetograms. The necessary 
conditions of the force-free field are a lso expressed as: \F X \/F p <sc 1, \F y \/F p <§: 1, and \F Z \/F p «: 1 



(Met calf et al. 



1995 



Moon et al. 



20021) . that is if three parameters F x /F p , F y /F p and F z /F p are so 



small that they are n egligible then the magnetic field can be regarded as a force-free completely. 



Metcalf et al. 



(I1995|) suggested that the magnetic field is a force-free if the FJF P is less than or 
equal to 0.1. The calibration coefficients Cl and Cj may affect the three parameters of F x /F p , 
Fy/F p and F-JF Thus, two sets of calibration coefficients mentioned above are applied to this 



study . Casel: C L and C T are chosen as 8381 G and 6790 G (|Su & Zhang 



2004 



19961) . respectively. Casell: Cl and Cj are chosen as 10000 G and 9730 G (|Wang et al. 



Wan g. Ai & Deng 



1996), 



respectively. 

Fig Q] shows the possibility density function (PDF) and scatter diagrams of F x /F p , F y /F p and 
F z /F p for the selected magnetograms (Casel). The mean values of absolute F x /F p , F y /F p and 
FJF P for all selected 925 magnetograms are 0.077, 0.109 and 0.302, respectively. The amplitudes 
of F x / F p and F v /F„ are evidently smaller than that of F z /F p , which is the same as previous 



study (|Metcalf et al. 



1995D . Therefore F z /F p can work as a criterion indicating a force-free or 



non-force-free field more evidently. From the distributions of PDF and scatter diagrams, the most 
of magnetograms have the amplitudes of F-JF p distributed outside the zone consisting the width 
of + 0.1, which means the most of the photospheric magnetic fields deviate from a force-free 
field. There are about 17% of the magnetograms with the values of F z /F p less than 0.1 for Casel 
(38% of the magnetograms with the values of F z /F p less than 0.2). To see the relation between 
F x /F p , Fy/F p and F z /F p and magnetic field strength, Fig|2]shows F x /F p , F y /F p and F z /F p vs the 
magnetic components of B x , B y and B z for the selected magnetograms (Casel), where B x , B y and 
B z are the average of the absolute values of all pixels for each magnetogram. It can be seen only at 
the bottom right of Fig|2]that the amplitudes of F z /F p decrease roughly as B z increase, and there 
is no evident correlations exist between the parameters (F x /F p , F y /F p or F Z JF P ) and magnetic 
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components, but according to Equation (flOT ). the value of FJF P should decrease as B z increasing. 

To study the effect of calibration coefficients on these three parameters of F x /F p , F y /F p 
and F z /F p , PDF and scatter diagrams of F x /F p , F y /F p and F z /F p of the selected magnetograms 
(Casell) are plotted in Fig [51 the amplitudes of F x /F p and F y /F p same as Casel are smaller than 
that of F z /F p evidently. For Casell, The mean values of absolute F x /F p , F y /F p and FJF p for 
all selected magnetograms are 0.078, 0.1 1 1 and 0.251, respectively, which are smaller than those 
of Casel on the whole. And also for Casell there are about 25% of the magnetograms with the 
value of F z /F p less than 0.1 (49% of the magnetograms with the values of F z /F p less than 0.2). 
It should be noted that even more F z /F p in the zone mentioned above (consisting the width of 
+ 0.1), most of magnetograms can not be regarded as a force-free. Besides, the widths of PDF are 
not narrow and the scatter of diagrams are diverge as well. In general, there is deviation of F z /F p 
from the zone. For Casell, F x /F p , F y /F p , F z /F p and magnetic components of vs B x , B y and B, 
are also plotted in Fig [4] to study their relations with the magnetic field strength. The results is 
consistent with that of Casel, only the the amplitudes of F z /F p decrease as B z increasing. 

To see the differences of results of F x /F p , F y /F p and F z /F p between two cases, the scatter 
plots of F x /F p , F y /F p and FJF P of Casel vs the corresponding ones from Casell are shown 
in Fig [5] The correlations are calculated, which are 0.994, 0.995 and 0.980 for F x /F p , F y /F p 
and FJF p , respectively. Also the linear fits (y = Ax + B) are done based on these scatter plots, 
the values of A are 0.0008, 0.0003 and 0.2034 and B 1.10, 1.10 and 0.95 for F x /F p , F y /F p and 
F z /F p , respectively. The correlation show there exists high consistence between these two cases. 
Nevertheless correlation of F z /F p between two cases is not so good as those of F x /F p and F y /F p . 
This may imply that the more attentions should be focused on the amplitude of F z /F p in order to 
understand the extent of force-free. 
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4. Discussions and Conclusions 

It is worth to study the force-free extent of the photospheric magnetic field while together 
with the case of magnetic field extrapolation, since it has been assumed that the coronal magnetic 
field is force-free and the photospheric one should be matched observationally. In this paper 
the results of force-free extent of the photospheric magnetic field are given through a statistical 
research, using 925 magnetograms corresponding to 925 active regions observed by SMFT at 
HSOS. 

A part of efforts to avoid the uncertainty of the data calibration is the employment of two 
sets of calibration coefficients to describe the force-free extent of the photospheric magnetic field. 
For the easel, the calibration coefficients C L and C T are 8381 G and 6790 G, the mean values 
of absolute F x /F p , F y /F p and FJF p for all selected magnetograms are 0.077, 0.109 and 0.302, 
respectively, and for the casell C L and C T are 10000 G and 9730 G, the mean values of absolute 
FJF P , Fy/F p and FJF p are 0.078, 0.111 and 0.251, respectively. There are 17% and 25% 
magnetograms with their value of F z /F p less than 0.1 for easel and casell, respectively, in other 
word, 17% of easel and 25% of casell are force-free. Although there are differences between two 
cases, between them the correlation are as high as 0.994, 0.995 and 0.980, for F x /F p , F y /F p and 
F z /F p , respectively. Consequently we concluded that large part of the photospheric magnetic 
fields do not belong to a force-free. So before extrapolating a magnetic field, the force-free extent 
of the photospheric magnetic field should be adequately considered. We note that F z /F p decrease 
with the increase of B z , which is a more important parameter indicating whether a magnetic field 
is force-free or not, since the amplitude of F z /F p is larger than those of F x /F p and F y /F p . From 
Equation (flOl) . it can be seen that F z /F p may be neglected when the amplitude of B z is comparable 
to ^B 2 X + Bj, while they are apparently different between B z and ^Jb x + B 2 , then the amplitude 
of FJF P should be enlarged. According to the Equations (|8l- ([T0T ) and observatories as well, the 
correlations between F x /F p , F y /F p and magnetic field component are not evident. 
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The strengths of active region magnetic fields observed by SMFT at HSOS are determined 
through calibration, under the conditions of a weak-field approximation and the linear relations 
between the magnetic field and the Stokes parameters /, Q, U and V. Additionally, comparing 
the recent data obtained from space satellite, the data of SMFT at HSOS has lower resolution and 
more uncertainty of magnetic amplitudes. These disadvantages may affect the statistical results we 
have acquired. However the statistical analysis associated with its results may have an significance 
as a practical reference, because SMFT has observed photospheric magnetic fields more than one 
solar cycles at HSOS and its reliability observations been studied strictly and adequately. It is 
hoped that the statistical results can be obtained based on high resolution data in the future. 
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Fig. 1. — PDF and scatter diagrams of F x /F p , F y /F p and F z /F p for the selected magnetograms. 
The mean values of absolute F x /F p , F y /F p and F z /F p are plotted and indicated by < \F X /F p \ >,< 
\F y /F p \ >, < \FJF P \ >, respectively. (Casel) 
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Fig. 3. — PDF and scatter diagrams of F x /F p , F y /F p and F z /F p for the selected magnetograms. 
The mean values of absolute F x /F p , F y /F p and F z /F p are plotted and indicated by < \F X /F p \ >,< 
\F y /F p \ >, < \FJF P \ >, respectively. (Casell) 
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Fig. 4. — F x /F p , F y /F p and F z /F p vs B x , B y and B z for the selected magnetograms. (Casell) 
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Fig. 5. — Scatter diagrams of F x /F p , F y /F p and F z /F p deduced from Casel and Casell, respec- 
tively, for the selected magnetograms. 



